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ABSTRACT 



We study the extent to which the cosmological fine-tuning problem - why the relic density of 
neutralino cold dark matter particles x is similar to that of baryons - is related to the fine-tuning 
aspect of the gauge hierarchy problem - how one arranges that Mw <C Mp without unnatural 
choices of MSSM parameters. Working in the minimal supergravity framework with universal 
soft supersymmetry-breaking parameters as inputs, we find that the hierarchical fine tuning is 
minimized for Q x h 2 ~ 0.1. Conversely, imposing VL x h 2 < 1 does not require small hierarchical 
fine tuning, but the exceptions to this rule are rather special, with parameters chosen such that 
m x ~ Mz/2 or M^/2, or else m x > m t . In the first two cases, neutralino annihilation receives 
a large contribution from a direct-channel pole, whereas in the third case it is enhanced by the 
large top Yukawa coupling. 
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One of the important philosophical issues to be addressed by any dark matter candidate is 
the extent to which a relic density of interest to astrophysicists and cosmologists, 0.01 < Qh 2 < 
1, is natural. Indeed, this question is often posed to proponents of favoured candidates such 
as a neutrino, the axion and the lightest supersymmetric particle, assumed to be the lightest 
neutralino x 0- I n t ne case °f a neutrino, the see-saw mechanism explains in a natural way 
why m„ <C m q or me, but does not lead inexorably to Q u in the interesting range. In the 
case of the axion, experimental and astrophysical constraints restrict its relic density to the 
interesting range, but this is not yet predicted by any deeper theoretical argument, although 
such an argument may yet be found 0. In the case of supersymmetry on the other hand, it 
is well known and will be emphasized below that there are indeed good physical reasons for 
expecting cosmologically significant relic densities ||. 

The essential motivation for supersymmetric particles to appear below the TeV scale is 
provided by the gauge hierarchy problem. Supersymmetry by itself does not explain why 
Mw <C Mp, but it does enable such a hierarchy to be stabilized against the effects of radiative 
corrections, averting the need for fine tuning and rendering the gauge hierarchy technically 
natural [|J. The appearance of supersymmetric particles at the TeV scale is also supported 
||, P by the experimental value of sin 2 9w, in accord with supersymmetric grand unified theories, 
and by the indications that the Higgs boson may weigh around 100 GeV [0 , in agreement with 
supersymmetric model calculations if sparticles weigh < 1 TeV ||. 

The TeV scale also arises as a possible characteristic mass scale for a cold dark matter can- 
didate. Particles that annihilate via conventional point-like interactions generically have Q ~ 1 
if their masses m ~ \/Mp x Tcmbr, where the cosmic microwave background temperature 
Tcmbr = 2.73 K @. It so happens that y/Mp x Tcmbr ~ 1 TeV, making it plausible that any 
relic with mass around the electroweak scale might have a cosmological density of astrophysical 
interest. 



In the case of the supersymmetric relic Xi detailed calculations have been performed fTOfl . 



and it has often been observed that 0.01<fi x <lisa generic feature of parameter choices in 
the minimal supersymmetric extension of the Standard Model (MSSM). Moreover, it has also 
often been argued that restricting fl x h 2 < 1 suggests very strongly that m x is at most a few 
hundred GeV |L1]. However, it is known that there are rays in parameter space along which 
Q x h 2 may be kept small even though sparticle masses grow large. The purpose of this paper 
is to bring together and complete these observations, with the aim of clarifying the extent to 
which the supersymmetric resolution of the fine-tuning problem is related to the cosmological 
fine-tuning problem. 



The criterion for hierarchical fine tuning that we use is that championed in [12, [L3| , namely 



the logarithmic sensitivity A of Mz to variations in MSSM input parameters a,: 

a - i a i a ai dMz m 
A = max|A ai |: A ai = — — (1) 

The lower bounds imposed on A by experiments at LEP and elsewhere have been discussed 
previously by several authors. There have been extensive discussions how this price varies as 
a function of tan/3, the ratio of Higgs vacuum expectation values (vev's), and how the price 



1 



may be reduced if some underlying theory imposes relations between some MSSM parameters 
1l4| , p~5| , Here we extend these discussions to include a cosmological dimension. 

We study the extent to which the cosmological fine-tuning problem is related to the fine 
tuning of the gauge hierarchy, as measured by the quantity Ao (0). Making an extensive search 
of the MSSM parameter space, we find that minimizing A leads to Q x h 2 ~ 0.1. Moreover, 
we find a general correlation between those parameter choices with larger Q x h 2 and those with 
larger A . On the other hand, imposing Q x h 2 < 1 does not exclude all models with large 
Ao- We can distinguish three regions in the neutralino mass m x (or, correspondingly, in the 
GUT scale gaugino mass parameter M1/2) in which Q x h 2 < 1 with large Ao are possible. For 
m x pa Mz/2 or m x pa (that is for M1/2 ~ 150 GeV) large pole-dominated s— channel 

annihilation cross sections give a low relic density in ways unrelated to the value of Ao- In the 
intermediate region, Mz/2, M^/2 < m x < mt neutralino annihilation proceeds mainly through 
slepton exchange, and £l x h 2 < 1 then implies m.o/Mi/2 ~ 0(1) and small A . For m x > m t , 
Q x h 2 < 1 is also possible in models with relatively light stop, so that the t— channel annihilation 
into ti pair is enhanced again leading to a lower relic density. Within the minimal SUGRA 
framework such models require large left-right mixing in the stop sector and hence, in addition 
to a large top Yukawa coupling, a large value of A . Thus, models with Q x h 2 < 1 and large 
fine-tuning are rather special: some have m x pa Mz/2 or M^/2 whilst others have m x > m t 
and relatively light stop. The former class of special solutions can be exhaustively explored 
by chargino searches at LEP 200, whilst the latter would be absent if \A \ < 1 TeV. Our 
study does not find that the cosmological and hierarchical fine-tuning problems are equivalent, 
but it does confirm that an interesting cosmological relic density is indeed to be expected in 
supersymmetric models that do not exhibit extreme fine tuning. 



Our study of these issues is based on the survey of MSSM parameter space made in jnj IB 
which we review briefly here. As usual, we denote the Higgs mixing parameter by /i, and 
we assume the conventional minimal parameterization of soft supersymmetry breaking in the 
MSSM, via a universal scalar mass parameter mo, a universal gaugino mass parameter M1/2, 
and a trilinear (bilinear) coupling A (Bo). We assume that /io, m , M1/2, A and B Q are 
the appropriate inputs at the GUT scale Mqut = 2 x 10 16 GeV which should be used 
in the hierarchical fine-tuning criterion (H). These parameters are renormalized down to the 
electroweak scale in the standard way, leading to an effective potential that breaks electroweak 
gauge symmetry spontaneously, with a calculable value of tan f3. 

Our procedure for surveying the MSSM parameter space is, for reasons of convenience, 
to choose low-energy parameter sets that respect the experimental constraints and yield an 
appropriate electroweak vacuum, as determined using the full one-loop effective potential in 
the MSSM ||17|| . Specifically, using the measured value of Mz, for each value of tan/3 and a 



given sgn(/J,), we scan low-energy values of the left-handed doublet squark mass tuq, the right- 
handed singlet up-squark mass mp and the CP-odd MSSM Higgs mass Ma that are allowed 
by the experimental constraints fI7|. We then use the renormalization-group equations to find 



the corresponding allowed values of the GUT input parameters. 

As experimental constraints, we take into account the precision electroweak data published 
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at the Moriond conference fllql , and require A\ < 4 in a global MSSM fit. We also in- 
corporate LEP lower limits on the masses of sparticles and Higgs bosons Another im- 
portant accelerator constraint is provided by the recently measured h — > 57 branching ratio 



2 x 1(T 4 < Br(B -> X s j) < 4.5 x 1(T 4 |f20|, which we treat as described in [|TJ]. We find that 
the GUT-scale parameters corresponding to successful choices with an upper cut of 1.2 TeV on 
the soft masses Ma and the squark masses, vary over the ranges \fio\ ^2 TeV, m ^5 1.7 TeV, 
M1/2 ~ 600 GeV, \A \ ^ 5 TeV and \B \ ^ 3 TeV. Note, in particular, that the upper bound 
on M1/2 (and, in consequence, on the chargino mass) visible on the plots follows simply from 
the cut imposed on the scanning procedure. 

Finally, we note that the range of the relic neutralino density favoured by astrophysics and 
cosmology is fl x h 2 ~ 0.1 (see, e.g., [plU), but we do not use this as a constraint in our analysis. 
Rather, our aim will be to explore the extent to which this is a natural outcome for successful 
MSSM parameter choices with small values of the hierarchical fine-tuning measure Ao (|l|). To 
this end, within the context of the supergravity-based MSSM, we calculate the relic density for 
each of the models considered in the fine-tuning analysis of [14], |I6|| . Thus, for every allowed 
set of the GUT-scale parameters obtained by our scanning procedure, we have a calculation 
of the relic density Q x h 2 in terms of all of the low-energy masses, which are determined from 
M1/2, ^0 and Aq for the same fixed values of tan f3 and sgn(n). 

We note that there are two recent refinements of the analysis of the MSSM parameter space 
and the dark matter density that have not been included in this survey. One is the latest 
implementation |22| of the requirement that our electroweak vacuum be stable against possible 
transitions to vacua that violate charge and/or color conservation. This requirement tends to 
exclude parameter choices with mo/Mi/2 < 1/2, which do not have exceptional values of either 
A or Q x h 2 . We have also omitted the possibility of coannihilation ]23| between x an d the tr, 
which is the next-to-lightest supersymmetric particle in a generic domain of parameter space. 
This is important when m x £ l.lm? R , which is the case only for a very small number of the 
parameter choices in our survey. When it is significant, it tends to enable points with larger 
hierarchical fine tuning Ao to have a cosmologically interesting value of Q x h 2 . 

Results for the case tan/? = 2.5 are shown in Fig. 1. The top left panel displays directly 
the correlation we find between the hierarchical fine-tuning parameter A and the neutralino 
relic density. Here and elsewhere, the eight-pointed stars represent parameter choices where no 
specific direct-channel annihilation mechanism is dominant. The five-pointed stars represent 
parameter choices where m x ~ M^/2, so that XX annihilation via the direct-channel Z° pole is 
dominant. Because the width of the Z° is relatively large, the effect of s-channel annihilation 
through Z°'s occurs only in a small patch of the parameter space. The open circles represent 
parameter choices where m x ~ M^/2, so that XX annihilation via the direct-channel pole of the 
lightest MSSM Higgs boson h° is dominant. In this case, because of the very small width for 
h°, the suppression of the relic density actually begins even when 2m x ~ 0.8M^ [|4[] and thus 



can cover a broader parameter volume. Finally, the dots represent parameter choices where 
m x ^ m t , so that XX ~^ & annihilation is important. 

We see in the top left panel of Fig. 1 that (i) the minimum value of Aq ~ 13 is attained for 
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Q x h 2 ~ 0.1, (ii) the minimum value of A increases gradually for larger values of Q x h 2 , (iii) 
apart from a few exceptional points, mainly with m x > m t , all choices with Q x h 2 < 1 have 
A < 100, and (iv) apart from choices where XX ~* tt dominates, when Q x h 2 > 1 there is a 
clear tendency for its value to be correlated with that of Ao- The origins of the exceptional 
choices with low Q x h 2 and large Ao are seen in the top right panel of Fig. 1. They have values 
of Mi/2 ~ 100 to 150 GeV, corresponding to m x ~ Mz/2 and/or Mh/2. We also see in this 
panel that the choices with annihilation into it correspond to M\/2 ^ 400 GeV, which are the 
largest in our sample. It is therefore not surprising that these correspond to some of the largest 
values of A , as we see in the top left panel. 



The remaining panels of Fig. 1 show how the value of Q x h 2 is correlated with the values 
of other MSSM parameters. In the middle left panel, we see a clear correlation between Q x h 2 
and the ratio mo/Mi/2, again apart from a few choices where XX ~ * Z° or h° or XX 
dominates. Apart from these choices, we see that Q x h 2 is minimized for mo/Mi/2 ~ 1- In the 
middle right panel, we see that f2 x /i 2 < 1 is consistent with relatively small values of \A \ for 
the m x ~ Mz/2 or Mh/2 points and the points with no dominant annihilation channel, whereas 
large 4 is required when m x > m t . This effect is even more pronounced for larger values of 
tan/? (see Fig. 2 and related comments below). 

It is worth emphasizing that the parameter choices dominated by XX ~~ * Z° annihilation 
will soon be explored directly by LEP, as will to a large extent those parameters which lead to 
annihilation via the light Higgs pole. These choices have m x ± ^ 100 GeV, and hence should 
reveal their secrets when the LEP center-of-mass energy is increased to ~ 200 GeV in the years 
1999 and 2000. 

We have carried out similar parameter studies for several larger values of tan/3 < 30, as 
well as for tan/? = 1.65, and now discuss some similarities and differences in these cases. As 
a general rule, the cases with tan/3 > 2.5 are qualitatively similar to the tan/5 = 2.5 case. 
In particular, it is always true that minimizing Ao favours Q x h 2 ~ 0.1, as seen in Fig. 2 for 
tan/3 = 10, for example. On the other hand there is no trend for the upper bound on Ao 
to be improved if one selects £l x h 2 < 1. We see again in the top left panel of Fig. 2 three 
distinctive sets of parameter choices which satisfy this condition. The one with XX ~^ Z° or h° 
annihilation, is unrelated to the values of A , as seen in the top left panel. The case with no 
dominant annihilation channel has small Ao- We also see there that the XX ~ * tt cases have 
distinctively larger values of Ao. For XX ~ > tt cases the correlation of Q x h 2 < 1 with large 
negative 4 or very large positive A Q is quite pronounced and, as explained earlier, reflects 
the necessity of a light stop, i.e., of a large left-right mixing. The asymmetry in the 4 values 
follows from the renormalization-group running of the mixing parameter down to low energies, 
where it is related to its GUT-scale counterpart through the equation [p5fl : 



A t = 4,(1 - y) - 0{2)M 1/2 (2) 

where y is the ratio of the top Yukawa coupling to its infra-red quasi-fixed point value. The 
correlation between Q x h 2 and the ratio m Q /Mi/2 is less pronounced, and m Q /Mi/2 > 1 is possible 
for Q x h 2 < 1 even away from the exceptional cases mentioned above. 

The situation is rather different for tan/3 = 1.65, as seen in Fig. 3. Here, we see in the 
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top left panel that the tendency to favour low Vt x h 2 is less marked, though Q x h 2 ~ 1 is still 
preferred. We again see that parameter choices with Q x h 2 < 1 normally have A < 200, with 
exceptions provided by choices with important direct-channel Z° and h° poles, as seen in the 
top right panel of Fig. 3. We recall that the five-pointed star choices with m x ~ Mz/2 will be 
explored exhaustively by the LEP runs at Eqm = 200 GeV. We also see in the panel f ) of Fig. 3 
that LEP Higgs searches will be able to verify or exclude this possible value of tan (3: it predicts 
that M/j ^ 96 GeV, whereas LEP 200 should have a reach extending beyond = 100 GeV. 

The results presented so far have been obtained with a scan over tuq, mjj and Ma up to 1.2 
TeV. We see that the highest values of M1/2 in this scan (and the corresponding values of m x ) 
are consistent with Q x h 2 < 1, but at the expense of increasing fine-tuning. It is interesting to 
extend the scan up to higher values of soft masses so that an absolute upper bound on Mi/ 2 
(m x ) is obtained from the requirement Q x h 2 < 1. This bound is shown in Fig. 4, where for 
tan/? = 10 we show £l x h 2 versus A and M1/2, scanning over rriQ, mu and Ma values up to 8 
TeV. In Fig. 5a we plot Ao as a function of m x for models which give Q x h 2 < 1. We observe 
that the bound for the heavier superpartner masses is weak, around 1 TeV, but is saturated 
only for large A . 

Let us now summarize the story so far. Minimizing the gauge hierarchy fine-tuning param- 
eter A favours values of Q x h 2 close to the range favoured by astrophysicists and cosmologists. 
Conversely, restricting VL x h 2 < 1 favours models with relatively low values of Ao, with certain 
well-understood exceptions, some of which may soon be probed by experiments at LEP. 

The question now arises how sensitive these observations are to the fine-tuning criterion we 



have used. In a recent paper fL6 |, we have studied the consequences of postulating some linear 
relation between a pair of the MSSM input parameters, and we now discuss their possible impli- 
cations for the cosmological fine-tuning problem. Fig. 6 displays the implications of assuming 
a linear relation between /io and M1/2, for the specific case tan/3 = 2.5. We see that the global 
minimum of Aj^ is significantly reduced, and that the preference for fl x h 2 ~ 0.1 is maintained 
and even enhanced, as compared with Fig. 1 where no parameter relation was assumed. We 
also see that the patterns of correlations between values of £l x h 2 and A for three different 
sets of parameters are maintained, modulo an approximate overall rescaling in the values of 
A . The picture if a linear relation between /i and A is assumed is somewhat different, as 
displayed in Fig. 7. The minimum of A^ is again reduced, as compared to the case with no 
parameter relation, but there is no longer any preference for fl x h 2 < 1: indeed, the favoured 
value is Q x h 2 > 10. Finally, returning to Figs. 5b and 5c we see, for tan/5 = 10, the values of 
A MfM and A^ for points satisfying Q x h 2 < 1 as functions of m x . The values A A / M < 10 are now 
compatible with m x < 500 GeV. 

We conclude by reiterating that there is a significant correlation between the amount of 
hierarchical fine tuning and the relic cold dark matter density in the MSSM. It is indeed 
"natural" that the supersymmetric relic particle have Q x h 2 ~ 0.1 to 1, which we consider to 
be an attractive feature of this dark matter candidate, as compared to massive neutrinos or 
the axion, whose densities have no obvious reason to fall within this favoured range. However, 
this cannot be regarded as a hard prediction of the MSSM. Moreover, fine tuning is always 
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a subjective argument, rather than a hard-and-fast mathematical argument, and one would 
immediately embrace any experimental discovery of even an "unnatural" dark matter particle. 
Nevertheless, we find this correlation between cosmological and hierarchical fine tuning an 
interesting supplementary argument in favour of supersymmetric cold dark matter. 
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Figure 1: The price of fine tuning and fi x h 2 for tan j3 = 2.5, as functions of various variables in 
the minimal supergravity model. The experimental constraints described in the text are included. 
The five-pointed stars (open circles) represent cases where m x ~ M z /2 (m x ~ M h /2). The dots 
represent parameter choices where m x <: m t . Eight-pointed stars represent parameter choices 
where no specific direct- channel annihilation mechanism is dominant. 




Figure 2: As in Fig. 1, but for tan/3 = 10. 
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Figure 3: As in Fig. 1, but for tan (3 = 1.65. 
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Figure 4: The price of fine tuning and Q x h 2 for tan/5 = 10, with the scan over mq, mjj and 
Ma extended up to 8 TeV. 
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Figure 5: The fine-tuning measures A 0; A« M and as functions of the lightest neutralino 
mass for tan/3 = 10. Only points satisfying Q x h 2 < 1 are shown. 
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Figure 7: As in Fig. 1 (tan (3 = 2.5), but instead assuming a linear correlation between A Q and 
jiQ. We show only points with < A . 
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